library(igraph)
##
## Attaching package: 'igraph'
## The following objects are masked from 'package:stats':
##
## decompose, spectrum
## The following object is masked from 'package:base':
##
## union
library(plotly)
## Loading required package: ggplot2
##
## Attaching package: 'plotly'
## The following object is masked from 'package:ggplot2':
##
## last_plot
## The following object is masked from 'package:igraph':
##
## groups
## The following object is masked from 'package:stats':
##
## filter
## The following object is masked from 'package:graphics':
##
## layout
library(ggnetwork)
df <- read.table("~/Downloads/higgs-reply_network.edgelist")
ig <- graph_from_data_frame(df)
degs <- data.frame(ind=(ig %>% degree(mode="in") %>% log(base=10) ))
degs$outd <- ig %>% degree(mode="out")
in_hist <- ggplot(degs, aes(x=ind)) + geom_histogram()
ggplotly(in_hist, text=ind)
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
## Warning: Removed 20332 rows containing non-finite values (`stat_bin()`).
ebs_df <- data.frame(abs=edge_betweenness(ig))
g <- ggplot(ebs_df) + geom_histogram(aes(x=log(abs,base=10)))
subg <- subgraph.edges(ig,1:500, delete.vertices = T)
tiny_g <- ggnetwork(subg, arrow.gap=0,layout=nicely())
tg <- ggplot(tiny_g, aes(x=x,y=y,xend=xend,yend=yend)) + geom_nodes()
tg <- tg + geom_edges(arrow=arrow(length=unit(4,"pt")))
tg
